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Abstract. Some basic questions about the liydrodynamical approach to relativistic heavy ion collisions 
are discussed aiming to clarify how far wc can go with such an approach to extract useful information 
on the properties and dynamics of the QCD matter created. We emphasize the importance of the coarse- 
graining scale required for the hydrodynamic modeling which determines the space-time resolution and the 
associated limitations of collective flow observables. We show that certain kinds of observables can indicate 
the degree of inhomogeneity of the initial condition under less stringent condition than the local thermal 
equilibrium subjected to the coarse-graining scale compatible to the scenario. 
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1 Introduction 

Hydrodynamic approach has shown to be very successful 
to describe the global and collective features of the pro- 
cess in relativistic heavy ion collisions. In particular, the 
behavior of elliptic flow parameter V2 as function of cen- 
trality and transverse momenta data is well reproduced by 
hydrodynamic models [1-5]. It should be mentioned that 
other approaches based on binary collisions of constituent 
particles like parton cascade in general yield smaller col- 
lective flow values than the observed values. This success 
of hydro approach leads to the expectation that, from a 
detailed hydrodynamic analysis of experimental data, we 
may be able to determine the thermal properties of QCD 
matter such as the equation of state (EoS) and transport 
coefficients, which can be then compared to those obtained 
from the lattice QCD (^QCD) calculations. In fact, many 
works have been and are being done in this direction. Of 
course, these are valid and important efforts that must be 
done to verify how far we can go with a given working 
hypothesis, that is, the validity of hydrodynamics. 

On the other hand, these successes brought us several 
new interesting questions and mysteries such as very early 
thcrmalization. The most crucial one is that why at all the 
hydrodynamic approach works so well for such a violent 
and almost microscopic; collisional process. It is commonly 
believed that the basic hypothesis of hydrodynamics is the 



validity of local thermodynamical equilibrium. That is, in 
every space-time points, all the thermostatic state vari- 
ables are well-defined and associate thermodynamic cqiia- 
tions are valid locally. If this is true, and the hydrodynamic 
analysis is giving a unique scenario, then we are led to con- 
clude that the thcrmalization time and correlation length 
are really extremely small for the QCD matter produced 
in the nuclear collisions. This is an important consequence 
for further understanding of the QCD matter at extreme 
conditions. 

Note that in the above statement, there are two impor- 
tant 'if's, that is, 'if really the hydrodynamics necessarily 
implies the local thermal equilibrium' and 'if the present 
analysis's can be considered as a unique solution for the 
observed collective behavior of the relativistic heavy ion 
collisions'. In this paper, we would like to address these 
'if's, and discuss what will be the experimental measure to 
estimate the validity of the hydrodynamic statements. In 
particular, it is fundamental to discover observables which 
arc sensitive to the local dynamics of the system. Wc or- 
ganize the present paper as follows. In Sec. II, for the sake 
of later bookkeeping, we make a brief review of the struc- 
ture of the relativistic hydrodynamics keeping the above 
two if's in mind, in particular, from the view point of the 
variational principle. Wc then discuss how the; c;oncept of 
coarse-graining can be introduced in the derivation of the 
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hydrodynamic equations. It will be shown that the ideal 
hydrodynamic description can be considered as an effec- 
tive model in terms of course-grained field variables. In 
Sec. Ill, wc discuss the compatibility between the coarse- 
graining scale and observables. We argue that most of col- 
lective flow signals, if averaged over events, are insensitive 
to the detailed space-time structure of the dynamics, as 
far as the coarse-graining scale is compatible with these 
observables. The genuine local hydrodynamic signal with 
local thermodynamic equilibrium requires high resolution 
on the space-time evolution, and such an information can 
only be meaningful for each coUisional event and not for 
the averaged observables. In this work, using a controlled 
set of initial conditions with different granularity, we look 
for observables which reflect the event-by-event (EbE) dy- 
namical information and consequently are sensitive to the 
initial state inhomogeneity. Section IV is devoted to the 
discussion of our results and perspectives. 



2 Relativistic Hydrodynamics 



where T is the 3x3 stress tensor. 

In general, the direction of the energy flow vector does 
not necessarily coincide with that of the matter current, 
that is, not always oc u^. In such cases, we write 

(6) 



where we choose q'^n^ = so that n represents the density 
of the conserved quantity N measured in LL frame and 
<7^ is the diffusion current with respect to this frame. For 
the sake of later convenience, we introduce the following 
projection operators. 



(7) 



and 



: gf^^ _ uV, (8) 

which can be used to decompose vectors or tensors system- 
atically into the parallel and perpendicular components to 
the four - velocity u'^. For example, the above Eq.(6) is ex- 
pressed formally as 



(9) 



In this section, we review the derivation of relativistic hy- 
drodynamics (mainly for the ideal case) and discuss the 
meaning of local thermal equilibrium in the context of 
coarse-graining procedure. 

First, let us consider a (classical) matter carrying some 
conserved quantity N , to which four-current density n^{x) 
satisfies the continuity equation. 



(1) 



In this note, to avoid unnecessary complexity for those 
who are not used to the "general relativistic" notations, 
we consider only the case of Minkowsky coordinate, x = 
{a;°,x} with the metric giju^ = diag{l, — 1, — 1, — 1}. The 
energy- momentum tensor T'^''{x) of the system also con- 
serves, 



d.T^^ix) - 0. 



(2) 



In the case that T^^'^ix) is a 4 x 4 symmetric matrix, e.g. in 
the absence of electromagnetic fields or spin variables, we 
can diagonalize at any point x. From the physical require- 
ment, we know that there is only one time-like eigenvector 
u'', with positive eigenvalue e [6], 



(3) 



and identify 

These projection operators can be used also to decompose 
the energy-momentum tensor. Due to the symmetry of 
T'^'', in LL frame, we can write as 



T^" = eA 



A^AfT^p 



(11) 



where the last term corresponds essentially to the 3x3 
stress tensor T in the local rest frame. We can further 
decompose the tensor T into isotropic irreducible parts of 
the local rotational symmetry as 

(12) 



T = ltr(T)l-F [T - |tr(T)l] 



where 1 is the 3x3 identity matrix. This decomposition 
results in the original energy-momentum tensor as 

T""^ = eAi^^ - PA^l n^"", (13) 



where 



Q/3 



1 A 

3 ^-L 



We can identify this above time-like eigenvector, with proper 
normalization, 

= 1, (4) 

as the four-velocity for the Lorentz transformation from 
the observable system to the local rest frame of energy 
flow. Consequently, the eigenvalue e is the proper energy 
density, that is, the energy density measured in the rest 
frame of the energy flow at the space-time position x. The 
rest frame of the energy flow is called Landau-Lifshitz 
frame (LL). In LL frame, the energy- momentum tensor 
becomes: 

re 01 



(14) 

Note that at this point P is the dynamical pressure which 
is the sum of the hydrostatic pressure and bulk viscosity. 
And 

= A^^AfT^fi - \ AfT^pA^"^, (15) 

comes from the non-diagonal elements of the local stress 
tensor, ni^" is a symmetric tensor, satisfying 



0. 



(16) 



T 



(5) 



At this stage, the original 4x4 symmetric tensor T'^'^ 
(10 parameters) is expressed in terms of m^, e, P and 11 ^'^ 
(10 = 3+1+1-1-5). Together with the conserved current 
n^, we have to know the time evolution of these 14 vari- 
ables when their initial values are given. When the time- 
evolution equations form a closed system within these 14 
variables, then we have defined the system of partial differ- 
ential equations of (dissipative) relativistic hydrodynam- 
ics. 
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2.1 Ideal Fluid 

In some special physical situations, the total number of 
variables reduces drastically. Suppose that the 3 eigen- 
values of the local stress tensor are degenerated. Then T 
becomes isotropic and consequently the tensor iT'*'^ van- 
ishes. Furthermore, if there is no diffusion of the conserved 
current in LL frame, then q^^ also vanishes. In such a sit- 
uation, the total number of variables reduces to 6. Since 
the conservation laws, Eqs.(l) and (2) furnish 5 equations 
among them, we need only one equation to close the sys- 
tem. Usually we introduce the so-called Equation of State 
(EoS) which establishes a functional relation among the 
local quantities, s, P, n, as 



P = P(£,n), 



(17) 



which completes the system of ideal hydrodynamic evolu- 
tion. Note that in such a system, the dynamical pressure 
is identified as the hydrostatic pressure. In other words, 
there is no bulk viscosity. More explicitly, for the set of 
5 independent variables, e,n,u, Eqs.(l) and (2) can be 
written in the form of time evolution equations as 



dt'' 



7-e = -(e + P)5X, 



d 

^dt'' = 

d f e + j 
dt\ n 



= -d'P, 



(18) 

(19) 
(20) 



where {d/dt)e = d^e + v^diS. Roman letters were used to 
denote spatial only indexes {i = 1, 2, 3). This set of equa- 
tions can be solved together with Eq. (17) for a given ini- 
tial condition for {s,n,v). In the equations above, 7 = u° 
is the Lorentz factor, = /"f and n* = jn are, respec- 
tively, the (three) velocity and the density in the global 
observer's system. From now on, we use the convention 
that the symbol * is used to distinguish the value defined 
in the global observational system of the corresponding 
quantity without * defined in the local rest frame. 



2.2 Variational Principle 

We see that the hydrodynamic equation for an ideal fluid 

can be obtained from the continuity equations Eqs. (1) and 
(2) with additional assumptions such as local isotropy of 
the stress tensor, null diffusion of the conserved quantity in 
LL frame and the existence of the EoS. All these additional 
conditions can be satisfied if the thermodynamical equi- 
librium is satisfied locally. For this reason, it is commonly 
believed that the success of hydrodynamic description in 
relativistic heavy ion collisions indicates that local ther- 
modynamical equilibrium is attained in these processes. 

In fact, the set of hydrodynamic equations formally 
constitutes a local covariant classical field theory. If we 
consider that relativistic hydrodynamics as in a covariant 



theory, then the local thermal equilibrium should be at- 
tained in each space-time point. This is somewhat a phys- 
ically contradicting c;ondition, because to attain the ther- 
mal cquilibriiim, we need a large volume and time. Fur- 
thermore, even if we admit that the local thermal equilib- 
rium could be the true answer for the success of the hydro- 
dynamic approach in the relativistic heavy ion collisions, 
we don't know yet any quantitative measure to determine 
the precision of the hydro predictions for the properties of 
the matter in question. Here, to further clarify these ques- 
tions, it may be useful to derive the equations of motion 
for an ideal fluid dynamics from a different point of view. 

Consider a small volume element AV* of the fluid at 
the position x, whose velocity is v. The relativistic La- 
grangian of a particle with the rest mass m is given by 



L = — m/7 



(21) 



where 7 is the Lorentz factor of the particle. The four- 
velocity of the particle is given by 



7 

7V 



(22) 



as usual. 

Considering the fluid element as a particle, the intrinsic 

volume is AV = AV* /j and the (average) energy density 
£ in its rest frame is related to the rest mass of this piece 
of fluid as m = s-yAV*. The Lagrangian of the fluid is 
then written as 

L = -y eAV* (23) 

— 

where the summation is taken over every fluid element. In 
the limit of small volume element, we can write 



/ 



d^xe. 



(24) 



We suppose that in each fluid element, the proper energy 
density e can be expressed as function of the proper den- 
sity n of the conserved quantity as 



e = e{n). 

By definition, the current, = nu^^ conserves. 



0. 



(25) 



(26) 



Let us then apply the variational principle for the action, 



7- 



xs 



(27) 



together with the constraints, Eq.(26) and the normaliza- 
tion condition of the four-velocity, 



xM" = 1. 



(28) 



Introducing Lagrange multipliers for the constraints, we 
get 



S J d^x[e + X{d^n^')+^{i 



1)]=0, 



(29) 
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for the variations of (n, u'^, A, ^). It can be shown [7] that 
the variational procedure leads exactly to the conservation 
of the energy-momentum tensor, 



^^.T'''' = 0, 



(30) 



where T^"' = {s + P)u^'u'' - Pg^"' with 



dn 



■ e. 



Thus, the variational principle with the action Eq. (27), 

with the constraints, Eqs. (26) and (28) leads to the ideal 
hydrodynamic equations, Eqs. (18-20). Here, we showed 
for simplicity the case of only one conserved quantity, but 
it is straightforward to generalize when there are more 
conserved quantities, say, {rij, i = 1, . . . ,r}. We find that 
everything is the same, except for the definition of the 
"pressure", Eq. (31) to 



Hi- s. 

on 



(32) 



In the above we showed that the set of ideal hydrodynamic 
equations can be derived from a variational principle with 
the action Eq.(27). There, the most fundamental condi- 
tion is that the local energy density can be expressed as 
function of other extensive conserved quantities. 



£ = e(ni, . . . ,nr). 



2.3 Variation in Lagrangian Coordinates 



(33) 



For the sake of later discussion, we will repeat the deriva- 
tion of hydrodynamic equations using the Lagrangian co- 
ordinates. Let us introduce the space coordinate system 
{X} at the initial time to. We can write the fluid density 
distribution at this instant as 

no = no(X). 

The fluid element initially at the position X moves ac- 
cording to the equation of motion and occupies the other 
location, say x = r{t; X). Here, r{t; X) represents the tra- 
jectory of the fluid element initially located at X. Thus, 
the fluid configuration at the instant t in the observer's 
system is written as 



When riQ is a smooth function^, by a suitable variable 
transformation, we can consider Hq as constant without 
losing generality [8]. 

The action Eq. (27) can be written as 



(31) with the Lagrangian density 



-Je 



(38) 



(39) 



Note that this Lagrangian density docs not depend explic- 
itly on the field r(X,t) but only on dr/dt and V^r. The 
symbol d/dt means the time derivative in time for a fixed 
X and Vx stands for the gradient with respect to X. 
To calculate the Euler-Lagrange equation. 



d dC 



dC 



= 



dtd{dr/dt) \diWxr) 
we use the following properties of the Jacobian J. 



(40) 



d{Vxr') 



= Vx ■ (Vxr^ X Vxr") = 0, (41) 



and 



dJ „ a 

^ . Vy = J r. 



(42) 



In the above, the components {i, j, k) of the vector r should 
be taken to be cyclic. Thus, 



d{Vxr) 



k 



d 1 



-e(n) 



dJ 



dn n 9(Vxr) 



(43) 



and 

dC 
d{dr/dt) 



n*(x, t) = l d3XnS(X)53(x - r{t; X)), (34) ^^at we get 



de , ^(1/7) 

TL 

dn d{dr/dt) 
dr de ^ 

d,t dn 

e + P 

no u 

n 



which can be rewritten as 
n*{x,t) ■- 

where 

J = det 



d (e + P 



J 



dr 



nS(X), 



,axyx=r 

The local proper density n is then expressed as 



1 * 
n = — n 

7 



7J' 



(35) 



(36) 



(37) 



dt 



u 



rVrP 



(44) 



(45) 



which is exactly the relativistic Euler equation, Eq.(20). 
Here, as before, we defined P as 



dn 



n — e, 



given Eq.(25). 

^ See the later discussion when no is not a smooth function. 
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An important fact is that, the canonical momentum 
field of r{t, X) is 



dC 



e + P 



d{dr/dt) n 
so that the Hamiltonian of the system is given by 



H 



dt 



= j d^x [{e + P)v?+e\ 
= j d^'^T^^ 



- L 



still obtain the equations for average densities of conserved 
quantities. For example, let us consider a set of (very 
large number of) particles, instead of a truly continuum 

(46) medium, which is also quickly moving. In this case, the 
density njj is a sum of Dirac delta functions so that any 
smooth relation with the energy density becomes mean- 
ingless. However, in most cases we do not have very precise 
resolution neither in space nor in time. Then we may in- 
troduce an averaged density distribution by introducing 
the smoothing function W/j(x) instead of 5 function in 
Eq.(34) and also take a certain time average, to define the 
smoothed density distribution as 

(47) n*(x,t) = J dt' J d^Xn*o{X)Ur{t' -t)Wh{\x-r{t';X)\y 



which means that the Hamiltonian density in x coor- 
dinates is given by the (00) component of the energy- 
momentum tensor, T'"' = {e + P)u^u'^ — Pg^^ , as ex- 
pected. 

Furthermore, the total momentum P of the system is 
given by 



P = / rfXnS 



/ 



u 



d^-xi {s + P)7U 



(48) 



which is nothing but the integral of the space part of = 
(e + P)m°w/. Thus, the expression from the definition of 
the Hamiltonian density is: 



r°° = (£ + P)u2+£ 

can be rewritten as 

TO°uo + T^'ui - evP = 0. 



(49) 



(50) 



Note that this is exactly the 0-th component of the eigen- 
value equation Eq.(3). 



2.4 Coarse Graining and Hydrodynamic Modeling 

In the derivations above, we have assumed that the rela- 
tion in Eq. (25) (or the general case Eq. (33)) is valid lo- 
cally. This is somewhat a very severe condition, depending 
on the definition of the densities. For example, when we 
consider the hydrodynamical description of the relativistic 
heavy ion reactions, differently from the normal hydrody- 
namic applications for macroscopic fluids, the typical fluid 
element size can not be taken too much smaller than the 
whole system and also the time scale of the collective mo- 
tion can not be much larger than that of the microscopic 
one. In such case fluctuations and inhomogencities from 
the microscopic level are expected to be very large and 
relations like Eq. (25) becomes difficult to be well-defined. 

On the other hand, the derivation of relativistic hy- 
drodynamics from the variational procedure suggests that 
we may apply similar method for certain averaged densi- 
ties over some statistical ensemble of fluid elements and 



(51) 

where Ur{t) and Wh{x) are smoothing positive definite 
functions (in mathematics are also called moUifiers which 
are approximations to the identity), satisfying 



Ur{t) 
Ur{t)dt 



0, \t\ > T. 
1, 



limT 



/ 

, j ip{t)Ur{t)dt = (^(0), 



(52) 
(53) 

(54) 



/ 



Wh{v) ^ 0, |r| > h. 
Wh{v)d^v = 1, 



lim^^o j ip{r)Wh{T)d''r = (^(0), 

where tp is an arbitrary well-behaved function. A typical 
example of U and W is the Gaussian distribution. The 
positive parameters r and h are the scales of the time 
and space resolutions, respectively. We can also define the 
3-current, 

r(x,i)= j dt' j d^Xnl{X)^{t',X)Ur{t' -t)Wh(^^-v{t';X)\) 

We can see that these averaged density and current satisfy 
the continuity equation, 

^^Jdt'J d'XnUyi)fUr{t'-t)WH{\^-r{t';X)\) 
= - j dt' j d^Xnl{X)^Ur{t' -t)Wh[\^-v{t';X)\) 
= j dt' j d^Xnl{yi)Ur{t' -t)^Wh[\^-v{t';X.)\) 
= - j di' J d'Xn*o{X)UAt' -t)^-V,Wh{\^-r{t';X)\) 

Using this conserved smoothed current and density, we 
can derive the four-current, 
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and the proper density, 

The smoothed four- velocity field is then defined as 

1 .„ 



n 



Once we defined the (smoothed) velocity field above, we 

can consider the local rest frame at the space time point 
given by a; = (t, x) in terms of a Lorentz transformation 
defined by 

A{u)u = 








We can also introduce the smoothed energy-momentum 
tensor in analogous way. Let 

be the energy-momentum tensor associated to the matter 
at the space position x = r(i?, t) at time t. From this, we 
introduce the smoothed energy-momentum tensor, 

(55) 

As before, we can show that 

d;T^^ = J dt'J d^^'Ur{t'-t)Wf,{\^-^'\)d'^T^''{^',t') 



if 



From this smoothed energy- momentum tensor T'"', we can 
calculate also the smoothed proper energy-density as 



(56) 



The smoothed proper energy-density defined in this 
way is an average of the energy density viewed in the 
rest frame of the matter fiow. The average is taken over 
all contributions from the matter within the range of the 
coarse-graining scale in space-time. 

In terms of hydrodynamic modeling, we take these 
smoothed distributions as the dynamical variables that 
represent the system. For one coUisional event, where ev- 
ery microscopic state is specified, we can calculate these 
hydrodynainical variables. Inversely, however, for a given 
set of hydrodynamical variables there may exist many dif- 
ferent microscopic configurations that give rise to the same 
set of hydrodynamical variables. Let us call f2 the set of 
all coUisional events with microscopic configurations that 
produce a given profile of the four-current at the ini- 
tial time to- In other words, i7 is the ensemble of events 
that corresponds to an identical hydrodynamical profile. 
If we calculate £ at a given space-time point x for each 
of the events in il, they do not coincide in general. This 



is so even for j^{x,t), if t > to. However, if the coarse- 
graining size is large enough, the ensemble of microscopic 
matter which contribute to the hydrodynamic variables 
at the point x can be sufficiently large enough so that we 
expect that s and n will distribute sharply around their 
respective mean- values, say e and n, due to the central 
limit theorem. Furthermore, since e and n are the average 
energy and matter densities belonging to the same fiuid 
element specified by the coarse- graining volume, we ex- 
pect that they will have a strong correlation, in such a 
way that they are functionally related 



£ = s{n). 



(57) 



When we neglect the fluctuations in e and j^, we may 
think of a physical model in terms of the averaged matter 

current which defines n* and velocity field v as the dy- 
namical variable. For this, we may take the model action 
as 

(58) 



The variational principle for this action establishes the 
ideal hydrodynamic equation of motion for the coarse- 
grained variables. This model describes the time evolution 
of the system belonging to a given statistical ensemble Q, 
that is, the hydrodynamic initial condition. 

When the fiuctuations of hydrodynamical variables within 
the ensemble /? are not negligible, the above ideal fluid 
modeling fails and we have to extend the hydrodynam- 
ical variable as stochastic ones. Even in this case, it is 
possible to construct the hydrodynamical model with the 
"\Qcriational principle. In Ref. [8] , it is shown that the vari- 
ational principle for the action 

/ = /(£,n*,v) 

extended to fluctuating (stochastic) variables leads to the 
Navier-Stokes equation in the non-relativistic case. As ex- 
pected, the fluctuation is directly related to the dissipation 
coefficients. 

Therefore in the present vision, the viscosity is an un- 
known model parameter related directly to the coarse- 
graining scale. It may take different values than that ob- 
tained, for example, from lattice QCD calculations, ex- 
cept for the limiting case of thermal equilibrium. In fact, 
the simulations by ideal hydrodynamics with fluctuating 
initial conditions give rise to a similar behavior of pt de- 
pendence of V2 to that of the viscous hydrodynamics [9, 
10]. 

It is interesting to note that, although the coarse-graining 
procedure defined through Eqs. (51) and (55) is neither 
local nor covariant, the resulting equation from the varia- 
tional principle is written as if the theory is relativistically 
covariant^. 



^ This situation reminds us the similar situation of the fre- 
quently used mean-field approximation of an effective field the- 
ory. Even if the final form of the approximation keeps the sys- 
tem of equations as if they are relativistic covariant, physically 
it is not covariant. 
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3 Event by Event Hydrodynamics 

As shown above, the ideal fluid description of heavy ion re- 
action can be regarded as an effective model for the coarse- 
grained variables for the flow based on the assumption, 
Eq.(57) together with the model action, Eq.(58). In this 
vision, the success of ideal fluid modeling in relativistic 
heavy ion collisions is, in short, reduced to the adequacy 
of the approximations contained in Eqs. (57) and (58) for 
collective variables. These two equations allow for much 
wider microscopic configurations than those restricted by 
the "local thermal equilibrium" condition (as defined in the 
introduction), usually considered as the necessary condi- 
tion for the validity of a hydrodynamic description. 

As we lengthy described in the previous section, the 
definition of hydrodynamic variables requires the coarse- 
graining procedure. The size of the hidden ensemble Q 
depends on the coarse-graining size. The larger the size 
of coarse-graining is, the larger the ensemble Q becomes. 
For larger i?, the two conditions, Eqs. (57) and (58) have 
larger chance to be satisfied and the probability of the 
success of ideal fluid model increases. On the other hand, 
for larger coarse graining size, we loose the resolution 
in the space-time recognition. In terms of hydrodynamic 
variables, we will not be able to see the inhomogeneities 
in density profile that has smaller wavelength than the 
coarse-graining scale. This affects directly the class of ob- 
servables that the model can describe. The extreme ex- 
ample is when the coarse-graining size is larger than the 
system size and time evolution. Then the ensemble Q can 
be regarded as the statistical ensemble of the whole sys- 
tem, and the resultant system is a simple fireball model. 
Therefore, the thermal model for the particle ratio can be 
considered in this category. 

The coarse-graining size is thus intimately related to 
the class of observables and the validity of hydrodynamic 
description. For some observables which do not require a 
precise space-time resolution of the system, hydrodynam- 
ics description with a large coarse-graining scale maybe 
sufficient. In such cases, the real local thermal equilibrium 
is not necessary for the siiccess of the hydrodynamic de- 
scription. In addition, the experimental observables are 
usually averaged over collision events classified in terms 
of configurations rather loosely defined, such as centrality, 
event plane, etc. This may mask the important informa- 
tion contained in each hydrodynamic evolution. For ex- 
ample, the recent studies [11,12] of hydrodynamics with 
fluctuating initial condition on the elliptic flow V2 indi- 
cate that the linear relation of integrated V2 to the initial 
deformation parameter £2 are quite insensitive not only 
to the granularity in the initial condition but also to the 
transport coefficients, when averaged over many events. 
This means that the event averaged integrated v^, has a 
poor resolution for the hydrodynamic evolutions, so that 
the coarse-graining size of the respective hydrodynamics 
needs not to be small. 

In order to claim that the real local hydrodynamics is 
valid, we need to have observables that refiect the gen- 
uine hydrodynamic profile in evcint by event basis. For 
example, remnant of a sharp shock wave propagation, if 



exist, would be a good observable which tells the possible 
coarse- graining size of the collective flow, since a shock 
wave [13-15] is a genuine local hydrodynamic phenom- 
ena. The shock thickness shoiild not be larger than the 
coarse graining scale of the collective flow. Another exam- 
ple of genuine hydrodynamic signal may be the emergence 
of Kelvin- Helmholtz instability in longitudinal direction as 
suggested in [16]. 

In any case, it is extremely important to study the hy- 
drodynamic observables in event-by-event basis, as pro- 
posed a long time ago [17 20]. Recently, thanks to the 
high quality experimental data, information contained in 
event-by-event fluctuations in hydrodynamic analysis and 
their physical origin are being more actively investigated 
[21-27]. 

The key point is that when we apply the hydrodynamic 
modeling, we do not know a-priori the coarse-graining 
scale suitable for the fiow variables in the real scenario^. 
Formally we can apply the hydrodynamic description even 
to very inhomogeneous initial conditions. But this requires 
a correspondingly small coarse-graining scale so that the 
hydrodynamic description breaks down eventually. Thus, 
the question is to know how far we can infer the initial 
condition inhomogeneity within the hydrodynamic mod- 
eling without knowing the coarse-graining scale actually 
permitted in the process. 

In order to clarify the question of the validity of hydro- 
dynamic models and the coarse-graining scale, it is essen- 
tial to find out the set of observables which carry the in- 
formation of the dynamical effects coming from the initial 
inhomogeneities. For this purpose, we take the following 
strategy. First, prepare initial conditions with different de- 
gree of inhomogeneities in a controlled way and study the 
fiow dynamics from these initial conditions within an ideal 
fiuid modeling. That is, we assume that the hydrodynamic 
coarse graining scale is smaller than the inhomogeneity of 
the given initial condition. Then, look for the final state 
observables which reflect directly the scale of inhomogene- 
ity in the initial condition. Once such observables can be 
identifled, we may look for these in experimental data. If 
experimental data show such observables which charac- 
terize the scale of initial state granularities, we may con- 
clude that the coarse-graining size which still permits the 
fluid-dynamical description should be of the order of cor- 
responding granularity scale revealed in the data. Further- 
more, observation of such granularity observables indicate 
strongly the granularities in the initial condition, which 
should be useful clues for the initial state of dynamics cre- 
ated by QCD. 



3.1 Effects of granularity in the initial conditions 

As previously discussed, we are not able to know a priori 

the coarse-graining scale required for the hydrodynamic 
description of heavy ion collisions. On the other hand, the 

^ Usually we estimate this in terms of Knudsen number or 
mean-free path, but they are useful when the system is com- 
posed of gas. 
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Table 1. Sets of parameters used for our simulations. The 
quantities in the header are: the tube width, cr, the total trans- 
verse energy content, Et, the number of tubes for a central 
collision (fe — 0), N"^^^^ and the peak energy density of the 
tubes, e. 
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Fig. 1. Mean transverse momentum dependency on the emis- 
sion time for different granularities. 



scale is intimately related to the space-time resolution. 
Therefore, if the experimental data reveals a signal that in- 
dicate the granularity scale of the initial condition in terms 
of the hydrodynamic variables, we can infer that the corre- 
sponding resolution, and consequently the coarse-graining 
scale is compatible to this resolution. In this sense, it is 
very important to establish such a measure which reflects 
the granularity scale in the initial condition within the 
framework of hydrodynamic modeling. For this purpose 
we investigate how the granularity of the initial condition 
affects the flow observables. 

In order to analyze the granularity we propose a para- 
metrization of the initial conditions based on the Glasma 
picture. Each event is described as a set of overlapping 
longitudinal tubes where the total energy content, Et, and 
their width, cr, are taken as parameters [28, 11]. The tubes 
are distributed according to a wounded nucleons distribu- 
tion with the constants chosen to describe Au-Au colli- 
sions at \/ snn = 200GeV/A. 

In table 1 we show all the values that were used dur- 
ing our simulations. For each parameter set, we computed 
the 2D boost invariant hydrodynamic evolution of 2000 
events. The impact parameter range was chosen as 6 =0- 
12fm and centrality classes were computed according to 
the final thermal multiplicity of pions, iV^. For the sake 
of simplicity, we compare the results of the extreme cases 
denoted A, B, C, and D (with decreasing level of gran- 
ularity). The other simulations are used later when we 
systematically analyze the parameter dependency. 

The granularity of the initial conditions generated by 
this model is controlled by the two parameters that give 
the shape of the transverse Gaussian distribution of each 
tube: a and Et- More conveniently, one may also think 
of the peak energy density of the tube, e = Et/{2TTa^). 
Higher granularities are obtained by higher e and smaller 

(7. 

In these calculations, we observe that there is a correla- 
tion between the averaged mean-p^ of final particles, (pt), 
and the instant of their freezeout, tfq. Figure 1 shows that 



the high-pf particles tend to be emitted earlier than those 
of low Pt, except for the very early times. Furthermore, 
this feature is enhanced by the granularity. 

The result in fig.l can be physically understood as fol- 
lows. The rise of (pt) at early times is due to the initial 
acceleration of the fluid whose transverse velocity is set 
to zero at tq = Ifm. On the other hand, the decrease 
of (pt) in later times indicates that the bulk of fluid is 
subjected to smaller pressure gradients than the fluid ele- 
ments present in the periphery - which freezeout earlier. 
Therefore, within the hydrodynamic picture, one may ex- 
pect that the information about the space-time evolution 
of the system can be extracted from the pt dependence of 
flow observables. 



3.2 NexSPheRIO 

Following the idea in the previous section, we propose to 
calculate the flow observables for separate pt regions. We 
chose two Pt regions denoted as low for pt =0-0.5GeV/c 
and high for 0.5-5 GeV/c. In this section, we present cal- 
culations using the NexSPheRIOmodel [29,18,9]. 

As argued before, we expect the initial geometry to be 
more strongly reflected in early emitted particles - which 
tend to have higher pt - than in late emitted particles - 
lower Pt ■ To check this hypothesis in a full 3D model, we 
have computed the flow and eccentricity phases, respec- 
tively, _ 

sin(n(/)) 

1 " 



arctan 



cos{n(j)) 



and 



1 



rdrdcf) e{r, (fyr^ sin(n(? 



iJ'n = — arctan 

^ ' rdrd(j)e{r,(f)r'^ cos {n(j>) 
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where n is the harmonic number, A is the pt range and 
e(r, (j)) is the energy density distribution of the fluid in the 
initial condition. We are interested in the relative phases 
between the low- and high-pt and ^n- For the sake of 
brevity, we will denote Atp^ = \tpf^ — tf'„|. 

In fig. 2 we show the distributions of Atp^^^^' (dashed) 
and ZiV2°^ (dotted). The top three plots show the distri- 
bution for 0-10%, 10-20% and 20-30% centrality bins, re- 
spectively. The distribution of Z\-02°* is peaked at Aip — 
but is very broad and spans over all the possible values 
0-7r/2. On the other hand, Aip'^^^^ is much more concen- 
trated around Aijj — 0. These results confirm our expec- 
tation base on the result of fig.l. We observe that the 
high-pt particles reflect better the properties of the initial 
condition when compared to the low-pt particles. 

The centrality dependency of the standard deviation of 



NeXSPheRIO Fluctuating IC [0-10%] 



the distributions, std{Z\V'' 



high(Iow) 



}, are shown in the bot- 



tom plot in fig. 2. Note that std{Z\-02"^^} (normalized by 
the maximum value of 7r/2) reaches a minimum of 0.33 for 
semi-central collisions (20-30%) and is consistently lower 
than std{Z\i/;^°*}. 

Now we use the fact that closely related to ^2 

in order to experimentally estimate — ^2] by sub- 

stituting ^2 by the experimentally measurable quantity 
^high^ The solid lines in the plots if fig. 2 represent the 



distribution Z\-02'^'^''°^ 



- U/,high 



low I 



Note that this 



quantity is in turn an experimental observable. 

As expected, the distribution of Z\-02'^'^''°^ (solid) fol- 
lows closely that of Aipl^"" (dotted). This indicates that 
one can extract experimentally the difference between the 
phase of the late emitted particles and the initial phase 
- by means of the early emitted particles. Furthermore, 
this difference depends on the details of space-time evolu- 
tion, in particular depends on the level of granularity of 
the initial condition. 



3.3 Granularity dependence 

In this section we perform a systematic study of the de- 
pendence of Zi-^j'^*^''"^ with respect to the granularity, by 
changing the model parameters a and Et- 

The distributions for the four parameter sets (labeled 
A, B, C and D in tablel, with decreasing level of granu- 
larity) are presented in fig. 3. For the most smooth case, 
D, there is no phase shift during the space-time evolution 
and almost all events fall in the central bin, A-ip = 0. On 
the other hand, for the case A, where the initial condition 
is the most granular among the cases we studied, there is a 
broad distribution of Aip^^^'^^^°'" . This result confirms that 
this observable is indeed very sensitive to the initial con- 
dition granularity. Furthermore, as in the NexSPheRIO 



A I high I low 

scenario, Aip2 



follows A^}]^^. 

, high|loi 



We also compute std{Zi'!/'2'^ ™} for all values in ta- 
ble 1. The result is shown in the first plot of fig. 4 as a 
function of the tube width, cr, for different Ef. There is 
a clear granularity dependence of the observable. Higher 
granularities - small a and large Et - produce broader 
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(explained in the text). 
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(HL), 



distributions of std{Z\?/'2"^^''™}- This suggests that the 
experimental analysis of this observable may provide im- 
portant information about the geometrical properties of 
the initial condition. 

It is also argued that higher harmonics would be sen- 
sitive to smaller angular variation, thus, within this ap- 
proach, we could expect a higher sensitivity to the gran- 
ularity. However, since the observable proposed here is 
also affected by the evolution dynamics and the different 
freeze-out times, the result might not show a clear de- 
pendence. The granularity dependence for the higher har- 
monics is shown in fig. 4. The higher harmonics have larger 
width than the elliptic case, showing that they are more 
sensitive for the space-time resolution, as expected. We 
also observe that, there is a difference in the dependence 
to the variation of tube energy content Et . In the case of 
odd harmonics, there is no difference between the results 
of different Et configurations, whereas, in the case of even 
harmonics, the curves are separate. This distinction be- 
tween the behavior observed for odd and even harmonics 
is still under investigation. 



4 Discussion and Outlook 

In this paper, we discussed the meaning of hydrodynamic 
modeling for the process of relativistic heavy ion collisions. 
We argued that it can be considered as an effective model 
defined by Eq.(58) for coarse grained hydrodynamic vari- 
ables. There, the validity of Eq.(57) is a key point when 
the average over the ensemble fl is taken. 

As mentioned, the most important issue is that for the 
hydrodynamic modeling, we don't know a-priori the valid 
coarse graining scale. The final form of relativistic hydro- 
dynamics after the coarse-graining is a set of equations 
of local fields and does not have any information of the 
coarse-graining scale. Even if the coarse-graining size for 
the validity of Eq. (57) should be large, we may model 
the system as if the coarse-graining size is very small. 
Furthermore, depending on the observable, they can be 
very insensitive to the fiuctuations in each event to ex- 
tract the detailed dynamical information [12]. Some ob- 
servables carry only information which are less affected by 
the coarse-graining size. On the contrary, if the behavior of 
an observable is well described by the hydrodynamic mod- 
eling, we can conclude that at least the coarse-graining size 
should be compatible with the observable. In this sense, 
the success of ideal hydrodynamics for the description of 
elliptic flow parameter tells us that the hydrodynamic 
coarse-graining scale must be less than that of specified 
by this observable. 

Some specific combinations of observables analyzed in 
event-by-event basis can carry more detailed information 
on the hydrodynamic evolution of the system and many 
efforts are now being done using higher harmonics [30, 
21]. In the previous section, we showed another example. 
There, if the hydrodynamic scenario is valid even for a 
very spiky initial condition, particles with high momen- 
tum can be associated to the early stage of the evolution. 
These particles are found to be emitted from hot spots 
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Fig. 4. Granularity dependence of std{Atp. 
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close to the surface of the system, since the time scale of 
expansion of such hot spot is very short. Therefore, we 
may consider the momentum of observed particles as a 
measure of time in hydrodynamic evolution. In fact, the 
correlations between event planes defined in momentum 
bins separated by the average pt seems very sensitive to 
the inhomogeneity scale of the initial condition, although, 
unfortunately the present analysis does not indicate still 
how to determine the inhomogeneity measure of the initial 
condition from the experimental data. However, if these 
observables in real experimental data, together with all 
other collective flow observables are consistent to a hy- 
drodynamic model, we could say that the coarse-graining 
size is at least smaller than the inhomogeneity measure 
determined in this way. 
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